Abstract. In the present paper, we show that eUiptical 
galaxies (Es) obey a scaling relation between potential en- 
ergy and mass. Since they are relaxed systems in a post 
violent-relaxation stage, they are quasi-equilibrium grav- 
itational systems and therefore they also have a quasi- 
constant specific entropy. Assuming that light traces mass, 
these two laws imply that in the space defined by the 
three Sersic law parameters (intensity Sq, scale a and 
shape v), elliptical galaxies are distributed on two inter- 
secting 2-manifolds: the Entropic Surface and the Energy- 
T— I \ Mass Surface. Using a sample of 132 galaxies belonging 
^^ ■ to three nearby clusters, we have verified that ellipti- 
^^ [ cals indeed follow these laws. This also implies that they 
' are distributed along the intersection line (the Energy- 
Q I Entropy line), thus they constitute a one-parameter fam- 
^\ ■ ily. These two physical laws (separately or combined), 
I allow to find the theoretical origin of several observed 
^^ ■ photometrical relations, such as the correlation between 
-^^ I absolute magnitude and effective surface brightness, and 
^ ■ the fact that ellipticals are located on a surface in the 
■TJ" I [logiJeff, — 2.51ogEo,logz^] space. The fact that elhptical 
C^^ ' galaxies are a one-parameter family has important impli- 
^1 i cations for cosmology and galaxy formation and evolution 
^sL . models. Moreover, the Energy-Entropy line could be used 
(<— ^ ' as a distance indicator. 

o- 

^^ I Key words: galaxies: clusters: individual (Coma, Abell 
O . 85, Abell 496) — galaxies: elliptical and lenticular, cD 
I ' — galaxies: fundamental parameters — distance scale — 
H ■ gravitation 
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1. Introduction 

Elliptical galaxies (Es) present a striking regularity in 
their global luminosity distributions in the sense that, 
within a wide range of sizes, their light profiles can be de- 
scribed by simple functions, such as the Sersic law (Caon, 



Capaccioli fc D'Onofrio 1993; Graham fc CoUess 1997; 



Prugniel & Simien 1997), which is a generalization of the 
de Vaucouleurs i?^/* profile. This regularity implies that 
ellipticals constitute a well defined family of galaxies. The 
regular properties of Es have been the subject of differ- 
ent approaches, concerning both photometric and spectro- 
scopic parameters, resulting in well known relations such 



as the Faber- Jackson and Kor mendy relations (Faber & 
Jackson 1976|; Kormendy 1977) or the Fundament al Plane 
( Djorgovski fc Davies 1987 ; Dressier et al. 1987 ). Es are 
supposed to be formed under non coUisional processes, 
where dissipation is expected to be negligible. Under these 
circunrstances, both the initial conditions and the gravi- 
tational forces are expected to have a crucial influence on 
their properties. 

In previous papers, we have addressed the question 
of how the properties of Es may be a consequence of the 



physical laws they should obey. In Gerbal et al. (1997) and 



Lima Neto, Gerbal fc Marquez (1999, hereafter Paper I), 
we have calculated the specific entropy of Es and shown 
that the results of analysing the cluster Es with Sersic 
profiles are compatible with a unique value for this specific 
entropy, assuming a constant mass/luminosity ratio. In a 
subsequent paper, Marquez et al. (2000| , hereafter Paper 
II) , we have shown that the specific entropy of Es was in 
fact correlated with the logarithm of the total luminosity 
while the specific entropy of simulated halos of coUisionless 
particles was correlated with the logarithm of the total 
mass. Moreover, we suggested that another physical law 



must be operating in order to explain that Es reside in 
a thin line embedded in a surface defined by a (almost) 
constant specific entropy. 

In this paper, we show that this second physical law 
is in fact a scaling relation between the potential energy 
and mass; photometrical relations can be naturally derived 
from this scaling law under some simplifying hypotheses. 
We also analyse what the consequences on the origin and 
evolution of Es are. In Sect. 2 we present a short the- 
oretical introduction to the concepts we will deal with, 
from which we derive a possible explanation for the ex- 
istence of an energy-mass relation, and we describe the 
passage from theoretical parameters to observational ones 
using the Sersic profile. We improve our analysis taking 
into account a possible dependence of the mass to light 
ratio on the total luminosity. In Sect. 3 we show that 
similar scaling relations are also found in other types of 
self-gravitating systems. In Sect. 4 we obtain the relations 
among the Sersic parameters and we compare the theoret- 
ical predictions with the results from observational data in 
Sect. 5. The discussion and conclusions are given in Sect. 
6. 

2. Theoretical background 

In the standard nrodel of structure formation, initially 
small perturbations in the cosmic density field start grow- 
ing in amplitude at the same rate. At a given epoch, when 
their overdensity relative to the average density is larger 
than a critical threshold that depends on the assumed 
cosmology, they stop expanding with the Universe and 
collapse. Bound objects then settle in an equilibrium con- 
figuration through gravitational and radiative processes. 



In the following (Sect. 2.1 and 2.2), we propose two 
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laws that elliptical galaxies should obey if they form and 
reach (quasi) equilibrium under the action of gravitational 
processes only. As a first approximation, we will assume 
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that radiative processes play a minor role compared to 
gravitation. 

2.1. Specific entropy 

Elliptical galaxies are thought to be in a quasi-equilibrium 
state, implying that they should obey the virial theorem. 
The second law of thermodynamics states that a system in 
equilibrium is in a maximum entropy configuration. Since 
elliptical galaxies are gravitational systems, they never re- 
ally reach an equilibrium state during the time scale of the 
two-body relaxation. However, even if the entropy 5' of an 
E galaxy is ever increasing on a secular time scale, after 
violent relaxation one may consider that the system is in a 
quasi-equilibrium stage, which is equivalent to saying that 
the entropy is quasi-stationary. 

Several works have been devoted to the problem of the 



entropy in gravitational systems (see for instance Mcrritt 
|(1999| ), and references therein). In a previous paper (Pa- 
per I), we have shown that, assuming that an elliptical 
galaxy is in a quasi-equilibrium stage, spherically sym- 
metric with an isotropic velocity distribution and constant 
mass-luminosity ratio, a thermodynamic entropy function 
can be defined. 

In Paper II, instead of using the thermodynamic 
definition of the entropy, we adopted the microscopic 
Boltzmann-Gibbs entropy, defined from the distribution 
function in the 6-dimensional phase space. In order to 
compare objects of different masses in a consistent way, 
we have introduced the specific entropy, i.e. the entropy 
normalized to the mass: 
_ S _ -Jflnfd^xd^v 

^~ M J fd^xd^v ' ^ ' 

where / is the distribution function. We have then shown 
that the specific entropy s, computed under the previous 
hypotheses and using the Sersic profile is almost the same 
among Es (within 10%). In fact, based on A^-body sim- 
ulations, we suggested that s is not strictly unique but, 
as a result of merging processes, weakly depends on the 
logarithm of the galactic mass. We then showed that the 
specific entropy and the luminosity are tightl y corr elated. 
Adopting an appropriate M/L ratio (see Sect. 2.3.1 below) 
this translates to: 

s = so + Ss\nM , (2) 

where sq and Ss are to be determined by observations. 

2.2. Scaling relation between potential energy and mass 

Let us consider a region with a characteristic length A^ and 
a mean density pi, in the initial density field (all quantities 
will be indexed with i - for initial). This region is supposed 
to be the location of a present E galaxy (its properties will 
be indexed with p - for present). 

One may suppose, at least to zeroth order, that the 
following quantities are conserved during the formation 
and subsequent evolution of the structure: 



1. the total mass: 

2. the total energy: 



M, = M„ 



M: 



Ei — Ep — E. 



The total energy at the initial time is: 
E^T, + U,^U^, (3) 

(where T and U are the kinetic and potential energies), 
because the kinetic energy is negligible compared to the 
potential energy. 

At the present stage, when the galaxy is already re- 
laxed, the two following relations are verified: 
Tp + Up = E, (4) 

and 

2Tp + Up^0. (5) 

Eq. (|^) is the virial relation, which of course does not hold 
in the i-case. 

Playing the game of addition, subtraction and com- 
parison, we finally obtain the usual relation for an ini- 
tially cold gravitational system that collapses and reaches 
equilibrium: 

Up ~ 2U, . (6) 

Notice that, if the initial kinetic energy is taken into ac- 
count, the above relation should read Up = 2Ui + 2Ti. 

In the initial linear regime all quantities, such as the 
typical perturbation length A.;, its mass Af, and its po- 
tential energy U, are related to one another by scaling 
relations, e.g. (we will adopt hereafter G = 1): 
M ex Af and U^ oc M'^/X, , (7) 

implying that: 

Using relation (0), we then obtain: 

Up oc Af ^/^ . (9) 

If we consider the (small) contribution of the initial ki- 
netic energy then the above relation will be tilted. Up ex 
7\,f '^+5/3 ^ "v^Jiere e is a small number depending on the exact 
initial conditions, e.g., the contribution of Ti to the total 
energy and the initial density fluctuation power spectrum. 

Equation (^ can be also written as: 
ln(C/p)-5/31n(M) = eo, (10) 

where cq is a constant to be determined by observations. 



2.3. From theory to observations 

2.3.1. Simplifying Hypotheses 

The relations given in Eqs. <M) and (00) are both theo- 
retically motivated laws that need to be translated into 
observable quantities if one wants to test them observa- 
tionally. In other words, we need to go from the observed 
two dimensional light distribution to a three dimensional 
mass distribution. Rather than doing a detailed dynam- 
ical and morphological model, we will assume some sim- 
plifying hypotheses: spherical symmetry, isotropic velocity 
dispersion tensor and a mass to light ratio independent of 
galactic radius R. 

The last of these hypotheses, i.e., AI/L = T{R) = 
constant, implies that the mass distribution follows the 
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light profile, with the same scale length and shape. Notice 
that the above "constant" may in fact depend on the to- 
tal galaxy mass or luminosity (as we will see below) . This 
hypothesis is clearly a rough simplification: X-ray observa- 



ness distribution of a sample of elliptical galaxies (see 
Sect, ^j), and to compute the physical quantities of inter- 
est, like the total mass, the potential energy, the specific 

in the Appendix for 



tion of giant ellipticals (e.g. Loewenstein & White 1999) 



gravitational lensing analyses (e.g. Grifhths et al. 199S), 
and dynamical studies (e.g. pcrhard et al. 2001 ), suggest 
that J^{R) increases with radius. This increase, however, 
seems to be important only beyond R ^ 2i?eff and current 
data is compatible with J-{R) ~ constant for R ^ 2i?off 
(iKronawitter et al. 2000|) . 



It has been argued (Berlin, Saglia & Stiavelli 1992) 



that a significant amount of dark matter can coexist with 
the stellar component without distorting the luminosity 
profile either if both the stellar and dark matter profiles 
have the same shape or if the dark matter is much more 
diffuse, i.e., with a larger scale length and lower central 
density. The former case corresponds to our hypothesis of 



entropy, etc. (see Tables |A.1| and | 
a summary of the formulae we use). 

Note that the Sersic law is a non-homologous general- 
ization of the de Vaucouleurs profile, in the sense that a 
third parameter describing the shape of the distribution 
(the structural parameter v) is introduced and left free, 
instead of being fixed to 0.25. This of course allows a bet- 
ter fit to the observations (Graham & CoUess 1997, and 
Papers I and II). 

We assume that only one dynamical component de- 
scribes an elliptical galaxy and ignore the possible pres- 
ence of nuclear cusps or embedded discs. This simplifica- 
tion is motivated by the fact that these components, if 
present, usually represent at most a few percent of the to- 



a cons ,ant T{R). In the latter case, the contribution of 



the da: k matter component will be important only in the 
external region of the galaxy. 

The spherical symmetry and isotropy hypotheses im- 
ply that the distribution function is simply a function of 
one variable, the binding energy. Although these hypothe- 
ses may hold in the internal region, they must be poor ap- 
proximations in the external regions of the most flattened 
galaxies. However, there is no clear correlation between 
ellipticity and luminosity and the effect of neglecting the 
flattening should not produce any systematic error. 

Finally, we will assume that mass and luminosity are 
related by a power-law: 

T{L) = kFL°', (11) 

where the index a may be zero (equal mass to light ratio 
for all ellipticals). The above functional form is suggested, 
e.g., by the analysis of the tilt of the Fundamental Plane 
assuming homology (e.g., Burstcin et al. 1997; l0rgensen 



fc Richstone 199 5 
et al. (1997| )~ 



tal galactic mass (c.f., ie Zeeuw & Franx 1991; Kormcndy 



and references therein; see also Saglia 



et al. |999| ), where a « 0.22-0.32. Notic e that the hier- 



archi cal mer ger simulations discussed by papelato et al 
(1995 |, |l997J) give a natural explanation for the origin of 
the Fundamental Plane under the hypothesis of a constant 
value of !F. 



2.3.2. Working model 

With the above hypotheses, the mass, potential energy 
and specific entropy are easy to compute from the light 
distribution. 



We stress here that we are doing a pure photometrical 
analysis and therefore such a working model is required 
(with all its simplifying hypotheses) so that we are able 
to derive dynamically related quantities like the specific 
entropy or the potential energy. 

3. Scaling relations: Elliptical galaxies only? 

The derivation of the scaling relation between potential 
energy and mass is based on fundamental principles (mass 
and energy conservation). It can therefore be applied to 
any object formed by the growth of density perturbations, 
driven by gravitational interaction. Dark matter (DM) ha- 
los are therefore expected to follow the same relation. 

Assuming that coUisionless A'^-body simulations pro- 
vide a suitable description of DM halo formation, rela- 
tion (|o|) should then also be found automatically in the 
results of such simulations. We have tested this hypothe- 
sis by considering a cosmological A'^-body simulation with 
512'^ particles in a 479 Mpc/h side box, characterized by 
the following parameters: a ACMD cosmological model, 
with density parameters f^o = 0.3 and Oa = 0.7, Hubble 
constant H — 100 ft. km s^^Mpc"^, h = 0.7, and normal- 
ization CTs — 0.9 (Yoshida et al. 2001). The correspond- 



Our working model is the Sersic (1968) luminosity pro- 
file: 

E = I]oexp[-(i?/a)'^], (12) 

which describes well the observed surface profile distri- 
bution of ellipticals (e.g., Caon, Capaccioli & D'Onofrio 



1993; Graham & CoUess 1997; Prugniel & Simien 1997, see 



also Sect. 6.1), and allows to compute analytically good 



approximations of all the needed quantities. We therefore 
adopt this profile to model the observed surface bright- 



ing particle mass is 6.8 x 10^^'h '''Mq. A sphere of ra- 
dius r — 7ft.~^Mpc around a very massive DM halo (M 
~ 2.3 X lO^^/i^i Mq, virial radius i?vir ^ 2.7/i-iMpc) has 
been selected in the original simulation, and re-simulated 
at higher resolution (a factor ^ 35 increase, to get a mass 
particle of about 2 x lO^h^^ Mq). Assuming a minimum 
mass of 10 particles per halo, this allows to resolve halos 
down to 2 X lOio/i-^Mg. 

Fig. M shows the distribution of potential energy ver- 
sus mass for the halos found within the selected sphere, 
their masses ranging from 10^^ to 10^^/i~^Mq. We have 
selected the halos that had their kinetic to potential en- 
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ergy ratio closer to the virial theorem value, i.e., the ones 
satisfying 0.8 < 2T/\U\ < 1.2. Due to the high resolution 
of the simulation, the computation of the potential energy 
is very accurate, at least for all halos more massive than 
~ 1O^^M0, i.e. composed by more than 100 particles. 




Mass (M^i^) 



Fig. 1. Potential energy-mass relation for about 29 000 
dark matter halos selected in a high-resolution A^-body 
simulation (see text). The solid line is the best fit to the 
data, with a slope of 1.69 ± 0.02, in very good agreement 
with the expected theoretical value of 5/3 (dashed line). 

A fit to the data indicates a power-law index of 
1.69 ± 0.02, very close to the theoretical value of 5/3. The 
difference may be due to the fact that the iV-body sim- 
ulation takes into account all dynamical processes (mass 
and energy exchange and/or loss) self-consistently. 

We have also verified that this relation is independent 
on cosmology. In fact, we found the same slope of the 
potential-ene rgy mass relation for D M halos in all the GIF 
Simulations ( Kauffmann ct al. 199£ ) , run for four different 
cosmological models: SCDM, OCDM, ACDM, rCDM. 

A c omparable simulation of dark m atter haloes was 
done by Jang-Condell fc Hcrnquist (2001 ) in a much lower 



mass range. In this simulation, each particle represents 
4.6 X 10** Mq, and the range of halo masses is 4 x 10^ Mq 
to 4 X 10^ Mq. The potential energy calculated for 2501 
bound halos also seems to follow a power law |C/| oc A/^/'^. 

This is an example of the universality of the scaling 
relation between potential energy and mass. Whereas out 
of the scope of this paper, numerous questions obviously 
arise, such as the dependence (or independence) of the 
index with redshift or with the cosmology assumed in the 
simulation. 

Interestingly, the analysis of ROSAT-PSPC images of 
24 clusters has also revealed the existence of a scaling law 
between the cluster potential energy and the mass of the 
X-ray gas, with a dependence U oc ]\,f^-T^±0-05^ again close 
to 5/3 (see Durret et al. 2001). Such a scaling law therefore 
appears to be general for self-gravitating systems. 



4. Predicted relations among the Sersic 
parameters 

We now show that the two relations introduced in 



Sects. 2.1 and 2.2, constraining the distribution of gravita- 
tional matter in galaxies, define two surfaces intersecting 
each other in the Sersic parameter space [p, a. Eg]. 

The first relation comes from the relation between the 
specific entropy, s, and the mass, M , that can be expressed 
as 

s^sa + Ss\nM ^ so + S^ \n{TL) , (13) 

cf. Eq. (H)q|. From the expression of the entropy in Ta- 
ble A.l, it is straightforward to obtain the Entropic Sur- 
face in terms of the Sersic parameters, Eq. (A.l) shown in 
the Appendix. 

The second relation is obtained using th e de finiti on of 
potential energy and mass given in Tables A.l and A. 2, 
that can be cast in a general form: 

U = kuM^ or, equivalently. In [/ = eq -f- /3 \n{TL) , (14) 
where cq = In kjj and we have dropped the index "p" . 
Using the luminosity and the potential energy derived for 
a mass distribution following the Sersic profile (cf. Ta- 
ble A.l), it is easy to obtain the Energy-Mass Surface in 
terms of the observed photometric Sersic parameters. The 
formula describing this surface is shown in the Appendix, 
Eq. ( [A^ ). 

A three dimensional representation of these two sur- 
faces is given in Fig. g, in the [v, a, Sq] space. The inter- 
section line (called the Entropy-Energy line in Paper II) 
is the locus on which elliptical galaxies are distributed in 
the space of the Sersic parameters. 



Entropic 
Surface 




In(Lo) 



.Energy-Mass 
Surface 



Fig. 

and 

[log So ) log a, i^]. In this example, we assume M oc L. 



2. 3-D representation of the specific entropy 
energy-mass 2-manifolds, using the coordinates: 



^ In Paper II, the specific entropy is defined with a different 
normalization; s' = — -jg- f dx'^dw^/ln/, where j dx'^Av'^ f = 
M. Here, we use the normalization j Ax'^dv'^ f* = 1 and com- 
pute the specific entropy a using the normalized distribu- 
tion function /*. With simple algebra it can be shown that 
s = s' +\nM. 
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Combining relations (A.l) and (A.2), i.e. calculating 
the intersection of the Entropic and Energy-Mass sur- 
faces, two-by-two relations between the Sersic parameters 
can be obtained. The general [a, ly] and [So,i^] relations 
are given in the Appendix. For the particular case where 
a = (constant M/L ratio), /3 = 5/3 (cf. Sectp^) and 
^s = (unique specific entropy), we have for the [a, i>] 

relation: 

3 13 

\n{a)+Faiiy) = --eo + -so + ^ InG, (15) 

withF,(;.) = ^lnR;{iy)-^lnM;{iy) + ^F2iiy), 
and, for the [Sq, v] relation: 
ln(Eo)+Fs(i^) = ^eo + iso-^hiG, (16) 



4 6 4 

withFs(i^) = ^-"gv 
The relation between a and So is obtained by combining 
relations (ll^) and (|l6|), with ^ as a parameter. 



-^lni?;(^) + ^lnM*(;.) + iF2(i.). 



5. Predictions faced to observations 

5.1. Data set and distribution of the v parameter 

The sample of elliptical galaxies we consider is described 
in detail in Paper II. It is composed of 68 galaxies in the 
Coma Cluster, 30 in AbeU 85 and 34 in Abell 496; their 
cluster membership is confirmed by their redshifts. 

For all these galaxies, we have fitted the observed sur- 
face brightness distribution (obtained by CCD imaging in 
the V band) by means of the Sersic integrated profile, i.e., 
using the growth curve fitting method (see Papers I and 
II, for a detailed description and discussion of the fitting 
procedure). As explained in our previous papers, SO galax- 
ies were excluded from our sample because they fall out of 
the correlations between the Sersic parameters. To be able 
to use all the data as a single set, the fitting parameters a 
and So have been expressed in physical units, i.e., in kpc 
and L0/kpc^ respectively^. 

While in the de Vaucouleurs law a single value of v is 
used {v = 1/4), our data fitting reveals that this param- 
eter can cover a range of values, which are not centered 
on 0.25. We give in Fig. the density distribution of v 



obtained by wavelet reconstruction [see Fadda, Slezak & 



Bijaoui (1998) for a description of the method]. The distri- 
bution is bi-modal, with a first maximum around v — 0.44 
and a second, less pronounced, one around v = 0.82 (close 
to 1, the typical value for dwarf spheroidal galaxies); no- 
tice that in the previously published i^-histogram (Paper 
I), this second maximum is already present, although not 
as significantly. This is due to the fact that we are now 
working with a larger sample, resulting from the combi- 
nation of three separate data sets. It is worth noting how- 
ever that our sample is more than 97% complete only up 
to magnitude 1/^ '^ 20 and the giant ellipticals of Coma 



were not included because of saturation of the CCD; be- 
sides, the CCD fields do not cover the entire clusters, but 
only the central parts (in particular for Abell 85 and 496). 
Therefore both the position of the two peaks and their 
relative intensities may be somewhat different when con- 
sidering complete samples. 




0.16 0.32 0.48 0.64 0.8 0.96 1. 12 1.28 1.44 
V 

Fig. 3. Histogram of the v distribution superimposed on 
the wavelet reconstruction of the density distribution. 



5.2. Entropy and energy-mass relations 



As discussed in Sect. 2.3, the mass, entropy and potential 



energy of all the galaxies in our sample were computed 
under the hypothesis that light traces mass. 

The uniqueness of the specific entropy of Es has been 
discussed and established in Paper I, at least at zeroth or- 
der. In Paper II we have shown that, as a result of energy 
and mass exchange in merging processes, sq weakly de- 
pends on the mass. A constant value for the mass-to-light 
ratio was assumed in these two papers. 

We now calculate the specific entropy as a function of 
mass for three mass to light relations: M/L oc i" (con- 
stant mass to light ratio), L"-^^ (from Fundamental Plane 
studies, e.g., Burstcin ct al. (1997 )) and L"'^ (extreme 



case). The results are shown in Fig. ^ and the best power 
law fits corresponding to Eqs. (|3|) and ( |l4| ) are given in 
Table 0. 

The potential energy is displayed in Fig. Has a function 
of total mass (adopting G = 1). A bisector-OLS (Ordinary 
Least Square) power law fit ( Feigelson & Babu 19921) gives 

FTl 



a indexes f3 = 1.79, 1.83 and 1.85 for M/L oc L",L°-^^ 
and L^-^ respectively (see Tabi c [l]) . However, if we impose 
(3 = 5/3 [as predicted in Sect. 2.2, see Eq. (^], the fit is 
still good enough (within 2a, cf. T abic p|) to c onfirm that 
Eq. (p^ is a good approximation. Fish (1964) analysed a 



We have used Ho = 100 km s^^Mpc"^ and go = 1/2 



set of 24 Es, fitting their surface brightness profiles with 
the de Vaucouleurs law, and found a power law index /? = 
3/2. The improved accuracy of both our data and our 
fitting procedure (using the Sersic profile), as well as our 
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10i» 10" 

Mass (solar units) 

Fig. 4. Specific entropy versus mass for three values of a. 
The bisector-OSL fits to our 132 galaxy sample are shown 
as full lines. 



larger sample, allow us to definitely exclude (3 = 3/2 (see 
also § §. 

Table 1. Fitting results for the specific entropy and po- 
tential energy scaling relations. The slope of the mass- 
luminosity ratio is a, sq and 6s are the parameters of the 
entropy-mass relation; eo and (3 are energy-mass param- 
eters. 



a 


So 


Ss 


eo 


/3 


0.0 


-8.0 ±1.2 


1.76 ±0.06 


3.68 ±0.03 


1.79 ±0.06 


0.25 


-6.5 ±1.0 


1.70 ±0.04 


2.74 ±0.05 


1.83 ±0.06 


0.50 


-5.5 ±0.8 


1.66 ±0.03 


2.08 ±0.06 


1.85 ±0.06 



Notice that the values of so for a = are different from the 
ones shown in Paper II due to our di ffere nt definition of the 
specific entropy (cf. Sect. W and Table A.l). 



5.3. Correlations between the Sersic profile parameters 

From a mathematical point of view the primary parame- 
ters a, So and v entering in the definition of the Sersic pro- 
file are independent from each other. But as discussed in 
Sect. PI, t he entropy and energy-mass relations, Eq. (A.l) 
and ( [A.2| ) , imply the existence of two- by- two relations be- 
tween the three Sersic parameters. From the observational 
point of view, we know that a, Sp and v are in fact corre- 
lated ( Young fc Currie 1995 ; Gerbal et al. 1997 ). We now 




10"' 101 

Mass (solar units) 

Fig. 5. Potential energy (with G = 1) versus mass for 
three values of a. The power-law fits (full lines) give slopes 
of 1.79, 1.83 and 1.85 respectively from top to bottom (see 
Table |l|. The theoretical relation, with an index of 5/3, is 
also shown (dashed lines). 



explain these correlations as originating from the physical 
laws discussed in Sect. |^. 

Figs, g and show the [a, ly], [So, H, and [a, Sq] rela- 
tions for our sample of galaxies. The corresponding the- 
oretical relations, computed for a=0.0, 0.25, 0.50 are su- 
perimposed on the data for two cases: the fixed value of 
/3=5/3 (the theoretical prediction as in Eq. g) and the 
values of P given in Table |[ The theoretical prediction 
lies close to the data points, in particular for big galaxies 
(those with the smallest v for /9=5/3 - see section |6.l[ ). 
We note that in both cases the effect of changing the M/L 
ratios is small. 



From theoretical to observed 



dations 



The well known photometric correlations which charac- 
terize elliptical galaxies, are usually expressed by means 
of astrophysical quantities like ReS, (A')eff ' and/or the ab- 
solute magnitude of the galaxies. Since these quantities 
are combinations of the primary parameters defining the 
Sersic law, if these observed photometric correlations do 
indeed contain information on the physics of the objects or 
on the processes that drove their formation and evolution, 
one should, in principle, be able to derive them from the 
"theoretical" relations, Eqs. (g) and (^0|). In doing so, we 
show that it is indeed possible to give a physical interpre- 
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Fig. 6. Pairwise correlation between the Sersic parameters [a, j/], [Sq, v] and [Sq, a] (counter clockwise from top left). 
The lines superimposed in each panel are the correlations predicted from the theoretical slope /3=5/3 for a^O.O, 0.25 
and 0.50. 
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Fig. 7. Pairwise correlation between the Sersic parameters [a, v], [Eq, i'] and [Eg, o] (counter clockwise from top left). 
The lines superimposed in each panel are the correlations predicted from the Entropic Surface (with a=0.0, 0.25 and 
0.50) and the fitted f3 values for the energy-mass scaling relations for these three values of a (see Table |l|). 



tation to some of the observed correlations and that they Before going through the comparison between theory 

are not just artefacts of the definitions of the parameters, and observation however, it is necessary to take into ac- 
count the following points: 
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The correlations proposed in the hterature (e.g., Young 
& Ourrie 1994; Binggeli & Jerjen 1998) resuh from 



fitting procedures; the coefficients are just empirical. 
Moreover, for the sake of simplicity, attempts are usu- 
ally done to linearize the correlations. On the contrary, 
the coefficients that we derive are the consequence of 
physical laws, and the relations we have found are non- 
linear. 
2. The intersection of the specific entropy and of the 
energy-mass surfaces is a line that we call the ^^Energy- 
Entropy line''' . The region of this line populated by 
ellipticals is not coplanar, but can be considered as 
located in a plane in a restricted range of values. 

5.4.1. Luminosity-effective radius relation 

From the intersection of the specific entropy and of the 
energy-mass surfaces, Eqs. (A.3) and (A.4) parametrized 
by r^, we can derive photometrical relations for Es galaxies 



using the formulae in Table A. 2. 

In Fig. @ we plot the total luminosity as a function of 
i?cff for our data. 

Using the definition of the potential energy, rela- 
tion (10) can be written as: 

eo = i[lnAf-31nrg], (17) 

where Vg is the gravitational radius and M the mass. As- 
suming that i?cff is proportional to Vg (which is a very 
good approximation for galaxies described by the Sersic 
profile, see Table A.2), we then have: 
logi = 31ogi?eff -f 31ogfc - logJ^+ 1.30288eo, (18) 

where T is the M / L ratio and k is the proportionality con- 
stant between i?off and rv. This relation is drawn with the 
dot-dashed line on Fig. g. Notice that the theoretical slope 
of the L-Rcs relation, in this case, depends only on the 
slope of the scaling relation between mass and potential 
energy. 

We can relax some simplifying hypotheses, proportion- 
ality between i?eft and r^ and constant M/L ratio and 
obtain a more general relation between L and RcS - the 
price to pay is a much more complicate expression, instead 
of Eq. (|8|) . The theoretical relations obtained for a fixed 
/3=5/3 (left) and for the fitted values of (3 (right) are su- 
perimposed for the three different values of a. Within the 
error bars, the relations are compatible with the data. The 
different curves vary with M/L when (3 is fixed (/3=5/3) 
but lead to similar curves when (3 is free. 

Our data and theoretical prediction are steeper than 



the results found in the literature. Binggeli et al. (1984) 
derive L oc R^g for galaxies with 0.2 <, Res ^ 5 kpc. For 
smaller galaxies, however, the slope steeps and reaches our 



theoretical value of 3. Schombert (1987) also finds that the 



L-Res relation is shallower for the brightest galaxies; for 
normal ellipticals he finds L oc Rig- Even though these 
low values are not incompatible with our data (within 
the error bars), it should be stressed that previous L-Rcg 



studies were based on the de Vaucouleurs profile, while we 
fit our galaxies with the Sersic law. As we will see below, 
different modelling (i.e., de Vaucouleurs or Sersic laws) 
of Es may change the observed correlations among their 
global parameters. 



5.4.2. Magnitude-mean surface brightness correlation 

Taking the definition of {fJ.)cS' 

{fi),s = -2.5 log L + 5 log i?cff + 2.5 log(2^) , (19) 

we can combine it with the energy-mass scaling relation, 
Eq. (|l7|), eliminating the effective radius. The results for 
the different values of a are plotted in Fig. |[ The theoret- 
ical relations obtained for a fixed /3=5/3 (left) and for the 
fitted values of /3 (right) are superimposed for the three 
different values of a. The different curves vary with M/L 
when (3 is fixed (/3=5/3) but lead to similar curves when 
P is free. 

If we make the same assumptions as in § p. 4.1 above, 
we obtain a relation between {fj.)cs and L, or equivalently, 
in terms of the absolute magnitude {M = — 2.51ogL): 

(A*)cff = iA^-51ogfc+^log.F-2.1715eo+^log(2^),(20) 

traced as a dot-dashed line in Fig. 0, where the variables 
have the same meaning as in Eq. ([l8|). Note that this re- 
lation does not predict two regimes but a single one. This 
is indeed in agreement with the observations when the 



hen 



data are analysed with the Sersic law, as done by ler 
Binggeli (1997), and as shown in Fig. g for our 



galaxies. Quantities like R^e or (/i)off are different when 



obtained through a de Vaucouleurs or a Sersic fit (Gra- 



bam & CoUess 1997; Jerjen & Binggeli 1997): this proba- 
bly explains why Binggeli et al. (1984), who used the de 
Vaucouleurs profile, obtain a broken power-law {fj,)oS~-M. 
relation, while Jerjen & Binggeli (1997) and us, using the 
Sersic profile, find a single power-law. 

5.4.3. Photometric plane 



Khosroshahi et al. (2000) have fitted a set of infra-red 
profiles of E and spiral bulges with the Sersic law. They 
have shown that these galaxies lie on a plane in the space 
defined by the set of coordinates [logi?eff, Mo,logi'], with 
fiQ — —2.5 log So- Notice that these authors use n — 1/v 
in the Sersic law. The equation they find for this plane is: 

(0.173 ± 0.025) logi?off - (0.069 ± 0.007)Aio = 

-log;/ -(1.18 ±0.05). (21) 

We will now derive a similar equation from a theo- 
retical point of view. For illustration purpose, we derive 
here the simplest case with constant mass to light ratio 
(a = 0) and (3 — 5/3. The general case is of course sim- 
ilar but more cumbersome. Including the definition of /io 
and the relation between i?off and a into Eqs. (A.l) and 
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fixed p = 5/3 
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Fig. 8. Luminosity versus effective radius for our 132 cluster ellipticals. The theoretical relations obtained for a fixed 
(3=5/3 (left) and for the fitted values of /3 (right) are superimposed for the three different values of a. The theoretical 
relation given by Eq. (Oq) is superimposed as a dot-dashed line. The continuous, dashed and dotted lines are obtained 
using the intersection of the specific entropy-mass and energy-mass relations, Eqs. (AT) and ( A. 2 ). 
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Fig. 9. Mean effective magnitude versus absolute magnitude. The results for the different values of a are plotted for 
the theoretical prediction /3=5/3 (left) and for the fitted values of (3 (right). The theoretical relation given by Eq. (20) 
is superimposed as a dot-dashed line. The continuous, dashed and dotted lines have the same meaning as in Fig. |8|. 



( [A.2| ), after some straightforward algebra, we obtain two 
equations: 

logii-cff - -7M0 - iogRlsi'^) + o loge (so - F^iiy)) , (22) 

15 9 

and 

logi?eff + ^Mo = logM*(^)-T941ogi?:ff(z.) 
5 

-1.303eo-T5113. (23) 

A linear combination between Eqs. (pa) and (£3) gives: 

0.1731ogi?eff- 0.069^*0 = -?^(i^), (24) 

with: 

K{iy) = -0.086125 logAf*(z/)-H 0.4262 log i?;:ff(i/) 

+0.1122eo + 0.025(so - -^2(1^)) + 0.13016 . 
Relation (24) is obtained from theoretical relations; on 
the other hand, relation (p^) is a fit to the observed data. 
The only differences between these two relations are the 
two right-hand sides, i.e. (— log;/) should be compared to 
K{v). Only the forms of the two functions are to be com- 
pared since the constant depends on the units chosen. In 
Fig. p^, we call attention to th e interval [0.25-0.7] in v, 
which khosroshahi et al. (2000| ) have used to define their 
photometrical plane. In particular, their relation given by 



Eq. ( |2l| ) and ours [Eq. (gj)] agree within 8% in this range. 
The good agreement between the observed points of Khos- 
roshahi et al. and our theoretical relation shows that the 
photometric plane may be understood as a consequence 
of the two laws discussed in the present paper. Note how- 
ever that it is possible to define a "plane" only because a 
limited range of values has been considered for v (see the 
beginning of § 5.4). 



6. Discussion and conclusions 

6.1. Ellipticals as a single parameter family 

The Sersic profile is certainly not the ultimate profile 
to reproduce the surface brightness of ellipticals, from 



brighter galaxies to dwarf ellipticals [see Jerjen & Binggel: 
|(1997D ]. The Sersic profile cannot apply to fine features 
such as boxy isophotes, central black hole, etc., but it al- 
lows a fit of sufficient quality for our aims, since we use 



it only to describe large scale properties (Caon, Capacci- 
Dli 



D'Onofrio 1993|; [Graham fc CoUess 1997| ; [Prugnie: 
fe Simicn 1997 ; Marquez et al. 200G|) . A finer description, 
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Fig. 10. The theoretical relation given by the generalized form of the right-hand si de of Eg. p4| (solid, dash ed and 
dotted lines) is superimposed on the fit (long-short dashed line) obtained by Khosroshahi et al. (2000) (right- 
hand side of Eq. ^l|) . The two vertical dotted lines correspond the range in v used by Khosroshahi et al. (200C| ) . The 
dash-dotted lines correspond to 8% errors on the fit. The shape parameter, n = 1/v, is given at the top of the figure 
to allow direct comparison with other authors. 



which is out of the scope of this paper, would certainly 
require another parameter to reproduce the whole profile. 

Although the use of the Sersic profile may seem es- 
sential to our calculations, it is nothing but an analytical 
expression which works well from the theoretical point of 
view (essentially because it provides a shape factor) and 
also observationally (it gives good fits). In fact any other 
function would be equally appropriate - provided it allows 
to calculate the entropy, energy, mass, etc. This is not the 
case for the de Vaucouleurs law which is characterized by 
two parameters only, and for this reason does not allow to 
define the specific entropy or energy-mass surfaces. Con- 
sequently, it is not possible to define the Energy -Entropy 
line using the de Vaucouleurs law. 

Notice that the de Vaucouleurs law is a homologous 
law, i.e., all models are deduced from a unique template by 
two scaling transformations. The introduction of a shape 
parameter in the Sersic law destroys - mathematically 
speaking - this homology. However, since all galaxies fitted 
by a Sersic law roughly have the same specific entropy, we 
can eliminate one parameter (Eg, say) by replacing it with 
by its expression as function of the two remaining param- 
eters, Eq = Eo(a, v). Therefore we have - mathematically 
speaking - only two independent parameters, like the de 
Vaucouleurs law. This procedure, however, does not lead 
to any kind of homology, for it is not possible to define 
a 'template' from which any Sersic distribution could be 
derived only by scaling transformations. 

Since all elliptical galaxies also obey the energy-mass 
relation, Eq. (pX)[), we can further eliminate another pa- 
rameter (a, say) with the result that a single parameter is 
sufficient to characterize a galaxy, for instance v. In this 
case every galaxy is defined by its shape parameter v. gen- 
eralizing the de Vaucouleurs law with the introduction of 
V allows to gather all ellipticals into a single large family, 
at least from the photometric point of view. 



The mass can then be rewritten formally as a function 
of V only as: 

Af(.)^^Eo(.)^rg), (25) 

where the equations for aiy) and Eo(i^) are given in 
Eqs. (na) and dlq). Moreover, very large galaxies are en- 
dowed with small values of v while dwarfs have large values 
of V (e.g., Caon, Capaccioli & D'Onofrio 1993, see also Pa- 
per I; the equation for a{i') together with the definition of 
RcS also implies this result). Thus, the parameter v plays 
the role of a concentration factor. This allows to under- 
stand why the de Vaucouleurs law cannot fit with good 
confidence these two extreme cases. 

Equation (^5|) above is in fact a bijection v -^^ M. 
During the processes of formation and relaxation, physics 
(described by the two laws discussed above) are necessarily 
acting so that at the end the above bijection is actually 
verified. 



6.2. Distance indicators 

The interesting possibility of using the correlation between 
[g, ly] as a distan ce indicator has been proposed by Young 
|fc Currie (1995|). This paper has been followed by con- 



troversies (Binggeli & Jerjen 1998; Young & Currie 1998) 



essentially due, as we understand it, to the following ques- 
tion: are distance indicators sufficiently reliable to be us- 
able for single galaxies, or can they only be applied to 
clusters globally? 

Since the correlations on which are based the indica- 



tors proposed by Young & Currie (1995) are consequences 
of the two laws discussed in this paper, it appears possi- 
ble to improve their quality. This is because the predicted 
relations are essentially non-linear, while the observed cor- 
relations have been tentatively fitted by phenomcnological 
linear laws, as we already pointed out. 
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This question has already been mentioned in Paper 
I and, in view of the predicted theoretical relations pre- 
sented here, we therefore hope to improve the distance 
indicator based on the shape profile of Es (to be analysed 
in depth in a forthcoming paper for a dozen clusters with 
rcdshifts between 0.06 and 0.37). 

6.3. Conclusions 

We have shown both from theoretical reasons and from 
observations, that elliptical galaxies obey a scaling rela- 
tion between the potential energy and mass (luminosity). 
In previous papers, we had already shown that Es share 
the same specific entropy (Papers I) with a logarithm de- 
pendence on the galactic mass (that may be due to an 
evolutionary process, e.g., hierarchical mergers Paper II). 

These two relations, Up-M and s-lnM, give an ex- 
planation to several observed correlations that have been 
proposed in the past by various authors, such as the cor- 
relations between the shape factor and a length scale, the 
correlation between the absolute magnitude and the cen- 
tral brightness, and the photometric plane. Therefore they 
constitute a theoretical background for a number of phys- 
ical properties of elliptical galaxies. 

The fact that elliptical galaxies lie on a line in the 
three dimensional space of the Sersic parameters implies 
that Es are indeed a one-parameter family. This has im- 
portant implications for cosmology and galaxy formation 
and evolution models. Furthermore, the Energy-Entropy 
line could be used as a distance indicator. 
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Table A.l. Various formulae related to the specific entropy and the energy-mass relation obtained for gravitational 
systems described by the Sersic law. The 2-D and 3-D refer to the dependence on Eg and po, respectively. The quantities 
M*, Mg and r* are defined in Table [A.2| . Note the different normalization used for defining the specific entropy from 
Papers I and II. 



Quantity 


Sersic expression 


Specific entropy 

2-D expression 

3-D expression 

Potential Energy 
2-D expression 

Specific entropy— Mass 
relation 

Energy— Mass 
relation 

2-D expression 

3-D expression 


s^ -j d-^rd^f* In/* ; j dVdS/* = 1 

s = 1 ln(.7^(L)Eo) + ^\na + F2(p) ; 

F2{v) ~ -0.795 \n{v) -^ + 3.85^/-''^ + InT [|] - 0.822 

s= l\n(T{L)jo)+G\na + Fi; 
\U\ =Ga^{2-KTT.ofxGl{u)- 


\i'r[(3-p)/iy]J ^ ■ ■ ' ' 
s = so + Ss\nM 

U = kuM^ 

eo = lnf/-/31nM 

3eo = In Eo - In a 4- In Ma (u) - 3 In r^ (u) + 3 In G 

3eo ^lnpo + In Mg* {v)-31nr; (u) + 3\nG 



Appendix A: Formulae 

General formulae useful for definitions or calculations of 
various quantities linked with our definition of the Sersic 



profile are given in Tables A.l and A. 2 



The relation between the specific entropy and the 
mass, Eq. (g) can expressed as a surface in the Sersic pa- 
rameters space with the help of the mass to light ratio, 
Eq. (|n|): 



a'S' In So = 5' In kp 



So 



+ 



a'S', + I 1 In 



27rr(2/i/) 



^;m + 



+ [2a'S'-- Una, 



(A.l) 



where a' = a + 1, S'^ = Sg — 3/2 and the function F2*(^) i^ 



defined in Table A. 2 



The scaling relation between the potential energy and 
the mass gives another surface in the Sersic parameters 
space: 

-a'(3'\nT,o = P' Inkp + \nku ~\n{4:7r^G) + 

+ [2{a'P' + 2)-3]lna, (A.2) 

where a' = a + 1, (3 ' = [3 — 2 and the function G2{i^) is 
defined in Table [O. 



The intersection of the above two surfaces yields the 
locus where elliptical galaxies are found in the Sersi c pa - 
rameter space. Eliminating Sq from Eq. (AJ) or Eq. (A.2) 
we get: 

{^ff + S',}\na = ,')'|,,„ - FJM] + 

+ 6'Jn[GATT^G;{iy)/ku]. (A.3) 

Note that In a versus v depends on the mass to light ratio 
only indirectly via 5a^ soj P and eg. Likewise, eliminating 
In a from the surface equations we obtain: 

(^/3' + j;)a'lnSo = (2a'/3' + 1)(F;H - so) 



[lp'-5'Aa' + l 



In 



2ttT{2Iv) 



- - 2a' 6' 
2 



^P' + 5'A\nki 



\n[G'iT:'^Gl{v)/ku] 



(A.4) 



The relation between a and Sq is obtained in paramet- 
ric form combining the above two equations and using v 
as a parameter. 
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Table A. 2. General formulae used in this paper. 2-D and 3-D quantities are marked with subscript "2" and "3", 
respectively. Quantities marked with a "*" are adimensional. Some of these expressions were first obtained in Paper I 
and are presented here for completeness. 



Quantity 


Theoretical Expression 


Analytical Approximation 


Sersic law 

2-D 


E(i?) = Eoexp(-(i?/a)'') 

V Vv \a J . 






3-D 




iM=io(i)"'cxp[-(r/a)n 
p ~ 1.0 - 0.60971/ -f 0.0546351/2 

L3(r)-ajo^7[ ^ ,(J J 




Mass profile 

2-D 

3-D 


M/L = T{L) = kpL" ; J^{R) = constant 
M{R) = T{L)L{R) 


M(r)=^(L)L3(r) 


Total Mass 

2-D 




■2- 






3-D 




M = fc^(a3joL^)"+^LJ = ^^ 


r3-pi 
// 




Total "magnitude" 


M = -2.5 log L 

= -2.5 log Eo - 5 log a -1- m* 


m* ~ -0.30413// - 1. 70786//- 1''''^'^^ 


Mean /i inside Rca 


(/i)cff = -2.51ogEo+m* + 

+51og(i?:ff)+2.51og2^ 




Radii 

Effcc. Radius 3-D 


Ri = a X R* 

7[(3-p)/j.,i?^,,ff] = ir[(3-p)/M] 


, „, 0.72701 -0.9877 In// „ „^„„, 
lni?3 eff ==; h 0.07021 


Effec. Radius 2-D 


7(2//., Rte) = \r(2/u) 


lni?:«. °-^0348 - 0.99625 ln//_^^g^22 


Gravit. Radius 


\U\=GM^/rg 


, . 0.82032 -0.92446 In// „„,,,, 
Inr: ~ -h 0.84543 

V 


rg versus _Reff 

Eo <^ jo 


a X Jo = Eo-^ 

., '' r(2///) 

""^"^^""2r([3-p]H 


Inr; ~ 1.16 -^ 0.98 In i?*H 

a jo — EoX 

X (0.076685 + 0.3253// - 0.041245//^) 



